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^■f^ , An efTective spin relaxation mechanism that leads to electron spin decoherence in a quantum dot 

' is proposed. In contrast to the common calculations of spin-flip transitions between the Kramers 

I doublets, we take into account a process of phonon-mediated fluctuation in the electron spin preces- 

■ sion and subsequent spin phase diffusion. Specifically, we consider modulations in the longitudinal 

' g-factor and hyperfine interaction induced by the phonon-assisted transitions between the lowest 

electronic states. Prominent differences in the temperature and magnetic field dependence between 
the proposed mechanisms and the spin-flip transitions are expected to facilitate its experimental 
verification. Numerical estimation demonstrates highly efficient spin relaxation in typical semicon- 
ductor quantum dots. 
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^ . Recently, there has been much interest in the study of electron spin dechoherence mechanisms in quantum dots 
I ' (QDs) since they are a natural candidate for the qubit operations of quantum computing. A typical approach to 
. this problem is to calculate the spin transition probability caused by the electron spin interaction with a thermal 
' bath [1-3] . A thermostat composed of nuclear spins is not an effective dissipative system due to the long nuclear spin 
^ . relaxation: however, the nuclei can assist in the phonon relaxation [4-7]. 

I Since spin-lattice relaxation (SLR) is assumed to be associated with spin-flip transitions between different electron 
. spin states, the non-diagonal spin-flip matrix elements must be taken into account. Most studies of SLR in the 
I literature have analyzed various mechanisms responsible for such transitions. Specifically, there are two approaches 
. to this problem. The first is based on the spin-orbital interaction which is particular to III-V [1] or Si/Ge [2,3] QDs. 
S ' The second approach incorporates the nuclear hyperfine interaction (HFI) [5-7] as a factor, leaving the time-reversal 
symmetry for the electron spin Hamiltonian [4]. 

Actually, there are no principal differences between the aforementioned approaches to the SLR in QDs and the early 
pioneering works performed in the 1960's for shallow donor relaxation (with the exception of the specific electronic 
energy structure and the infiuence of strain). Both the early works [8-10] and the later studies of QD SLR [1-3] deal 
. with the longitudinal (or energy) relaxation accompanying the exchange between the Zeeman and phonon reservoirs. 
' Note, however, that quantum computing is qualitatively limited by other relaxation processes which result in the 
. destruction of the electron spin phase coherence that can occur without energy relaxation. Development of transversal 
'if^ (or phase) relaxation was not the goal of the early studies nor the recent SLR research in QDs. 

. In this paper, we show that there are decoherence mechanisms specific to phase relaxation, which can be rather 
' more effective than energy relaxation. To make this assertion more clear, let us consider an electron spin s under 
"p^! . the influence of a magnetic field directed along the z axis with a randomly fiuctuating strength. In this case, the 
' projection of electron spin on the z axis Sz is conserved and no longitudinal relaxation occurs. Nevertheless, the phase 
^ , of electron spin will change randomly with the Zeeman frequency fluctuation Sfl resulting in a decoherence rate of 
I ^ (5f2^Tc. Here, Tc is the correlation time of these fluctuations [11]. This leads to a very important conclusion 

''2 ' that the phase relaxation time may not be determined by the matrix elements between the Kramers doublets, the 
Q ■ most important restriction on the longitudinal SLR. 

^ In order to calculate the fluctuations in the Zeeman frequency 17, we will consider the phonon- induced transitions 

between the lowest electronic states |fc) which possess different spin splitting and are not linked by the time reversal 
operator. The resulting electron spin phase diffusion due to the spin precession in a fluctuating fleld is expected 
to be an efficient relaxation mechanism if (i) the electron signiflcantly changes spin precession frequencies when the 
transitions between different states occur and (ii) the mean time Tc of phonon assisted transitions between these 
states is less than spin relaxation time T2 [11]. In the case of a QD, these two conditions can be easily realized due 
to shallow energy levels and the g-factor or hyperflne constant dependence on the orbital electronic states. 

We begin the quantitative analysis by deflning the Hamiltonian H over the basic functions reduced to a few lowest 
electronic states \k), which are involved due to the phonon-assisted transitions. We also assume that the single-electron 
problem in a QD without Zeeman energy and HFI gives the doubly degenerate energy spectrum Ek with eigenstates 
\k). It is conveniently (but not necessary) assumed that the spin splitting \n\ is small with respect to the energy 
intervals \Ek — Ek'\- With this assumption, the total Hamiltonian takes the form 

H — Hg + He + Hpfi + He-ph- (1) 
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The first term Hg = is the spin (or pseudospin) energy Hamiltonian; in the most general case, the frequency 
operator reads 

n = Y,^k,k'\k){k'\, (2) 

k,k' 

whore ^k,k' are the matrix elements of the effective field (in units of energy) taken between the \k) and \k') states. 
The spin-independent electron energies describe the Hamiltonian Hg 

H,=Y,Ek\k){k\. (3) 

k 

The Hamiltonians of the lattice and electron-phonon interactions have the usual form 

Hph = Y^i^q{a'laq + ^); (4) 
ife-pA= 13 S^,fc,|fc)(/c'|(a+ + a_,). (5) 

Here q = {q, h} represents the wave vector and polarization of a phonon with energy u)q] —q = {—q, x}, a+ and aq 
are the phonon creation and annihilation operators, j,, is the matrix element of the electron-phonon interaction, 
which depends on the material parameters and the geometry of the QD. In Eq. (1), the last three terms constitute 
the Hamiltonian of the dissipative sub-system responsible for electron spin relaxation, Hd = He + Hph + Hg-ph. 

We are interested in the evolution of electron spin s{t) = Tr{p{t)s\ (pit) is a density matrix) in a system with the 
Hamiltonian of Eq. (1). This evolution was shown to be described by the quantum kinetic equation (see Ref. [11]) 

j^s{t) =LJXs{t)-T {s{t) - so) , (6) 

where uj = ((l^ is an effective magnetic field, (...)= Trje^^''/-^ . . .}/Tre^^'^^'^ , and T is the temperature: uj and T 
are in units of energy. The matrix T of relaxation coefficients is composed of Fourier transformed correlation functions 
7^^ = 7^,^ {uj) = {dil^ (r) S^l,^)^=^ J ((5il^ (r) Sfli,) e'^'^dr; p, v = x, y, z; = il — /il). With a provision that the 

— oo ^ ' 

correlation functions are symmetrical, Jui, (w) = (w), the matrix T — tt ||r^^|| has a simple form in the frame of 
references related to the direction of the effective field uj \\ z: F^^ = 7°^ +njyy, Tyy = -y^^+njxx, ^'zz = ^'{-fxx + Jyy), 
F;, = -rn,u., (a* ^ v), where 7°, = 7,,(0), n = n{Lu) = (l + e^/^) /2, sq = -i tanh(c.V2r){0. 0, 1}. 

Thus, the problem of spin relaxation is reduced to the calculation of the correlation functions of the effective 
field operator with the Hamiltonian Hd of the dissipative subsystem. These calculations strongly depend on the 
specific form of Hd, the energy spectrum, and the quantity of electron states considered. Keeping this context in 
mind, we consider the more simple problem of electron fluctuations between only two discreet states |A;) = \g) or |e) 
corresponding to the ground state and the first exited (by interval 5) electronic energy level with Zeeman frequencies 
0.3 and Vi'^. Such a simplification allows us to easily perform all the necessary calculations in an analytical form. In 
addition, most of the important physics of the new mechanism under consideration can be obtained in the framework 
of this two-level model. 

Hereafter, it is convenient to introduce Pauh matrices ai, (72, 0-3 on the basis |e), \g), where according to definition, 
ci, (72, (73 are invariant with respect to the coordinate system rotation in contrast to actual spin matrices s. In terms 
of Pauli matrices, the Hamiltonian of dissipative subsystem takes the form Hd = Hph + das + ^qBqCTi (a+ -|- a_g). The 
electron spin Hamiltonian now assumes the form Hg = |(f^^ + 0^)5 + |(j1^ — 0^)ct3S, which defines the fluctuating 
part of the effective field 60, = ^{O^ — 03){a3 — ((73)) with (0-3) = — tanh((5/2T) and gives the correlation functions 
in Eq. (6) in the form 

7^. iu^) =\{%- %){^t - ^DUT); (7) 
UT) = {{a,{T)-{as)){a,-{a,)))^. (8) 

Skipping the derivation, we present the final result for the case w <C (5, 
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is the lifetime of the excited electron 



= 277^ (2n, + 1)5 K - 5) , (10) 

where Uq = {a^aq) is the phonon population factor for mode q. The parameter Tc has the simple physical meaning 
of the correlation time caused by phonon-assisted transitions between the g- and e-states. 

Actually, Eqs. (7)-(10) describe the problem under consideration in a very general form. Before we specify the 
electron spin interaction, which fluctuates due to the phonon induced transitions between \g) and |e) states, we 
provide the analysis of the SLR temperature dependence. In doing so, we note that the correlation time Tc given in 

Eq. (10) can be written as ts tanh {S/2T), where ts = 27rEq \Bqf S {ujq — 5) 

state with respect to the transition to the ground state through phonon emission in the limit T — > 0. Thus, at zero 
magnetic field (i.e., w = 0) the temperature dependence of SLR is reduced to 

Jo(T) = ^F^^^; ^'(a;) = (l-tanh^a;)tanha;. (11) 

When the magnetic field is strong enough, this factor also describes the temperature dependence of the transversal 
relaxation because only 7°^ survives in the matrix T [Eq. (6)] in the limit w^r^ ^ 1. 

As shown in Fig. 1, the pronounced maximum in the temperature dependence of F {S/2T) around T = S has a simple 
physical meaning. The left of the peak corresponds to the reduced hopping from the \g) to |e) state that decreases the 

difference — ^^^^ (or the amplitude of fluctuations). So, in the limit T <C 5 the fluctuations are frozen out and 

our mechanism becomes non-effective as exp(— 5/T). The negative slope in the right side (high temperature) arises 
due to the well-known effect of fluctuation dynamical averaging, which becomes more pronounced with an increase in 

temperature . 

As noted above, the correlation function tt^^^ = f (0| - f^f)2 Jo(T) describes the rate T''^ of the spin relaxation 
if longitudinal fluctuations dominate over transversal ones, 7^^ 3> j^xi^), lyy{^)- Let us apply the general theory 
discussed above to the mechanism of phase relaxation, which stems from the hopping between excited and ground 
states with different g factors. In the most general case, the reason for such a difference is the (7-factor dependence 
on the energy separation between the electron's discrete level and nearest spin-orbital split electronic band. For 
technologically significant III-V compounds, where the interaction with the valence band edge determines the deviation 
of electron (/-factor from free electron Lande factor go ~ 2 [12], one can find the amplitude of the fiuctuation Ag — 
5(.9o — g) (Ago + 2-Eg) / Eg (Ago + Eg), where Eg is the band gap, and Aso the spin-orbital splitting of the valence band; 
we also assume inequality 5 Eg. In the case of a silicon QD, Eg is a splitting i?i5 of A point in the Brillouin 
zone [9]. The relaxation mechanism due to the (/-factor anisotropy of the A band is not effective in the case of Si QD 
because of the specific valley orbital structure [3] . The final equation for the phonon-assisted rate of phase relaxation 
caused by the Zeeman energy fluctuations is given by 



i,P \E,(i\.„ + E,)) \2T 

where Hloq = g^sB, and /Ub is the Bohr magneton. One can see that our mechanism reveals a quadratic dependence 
of Tj"^ on the applied magnetic field B in contrast to the B^-dependence found in the previous calculations of 

longitudinal SLR [1 3,8,9] 

An estimation of excited state lifetime can be performed in terms of a deformation potential interaction and a 
model of lateral carrier confinement [3]. The matrix element of electron-phonon interaction between the \g) and |e) 
states in this model is Bq = iCy/hq/2pv^^VoJosc, where Jose = Josc{q) is a corresponding form-factor calculated in 
Ref. [3] , C is the deformation potential, and p, tiy and Vq are the density, longitudinal sound velocity, and volume of 
crystal, respectively. A straightforward calculation of inverse lifetime results in the expression 



'^5 ^ '30^2fe^,,2 



32 

where qs = 5/hv\\, a = 6/2meV?:, and me is a lateral effective mass. 



Jil- z'^)e-''^^-''Uz, (13) 
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To show the efficiency of the mechanism under consideration, we assume 5 = 2 meV [3] and calculate the relaxation 
parameters of Eqs. (13) and (12) for a GaAs QD under the magnetic field B = IT. We find that Tc — Tg = 5.8 x 10~^ 
s, and T2 = 2.0 s, 1.8 x 10~^ s and 5.4 x 10~* s for T = 1 K, 2 K, and 4 K, respectively. Similar calculations were 
provided for a Si QD with the same magnitudes of S, B and T: Tc = 3.5 x 10"'' s, and T2 = 2.6 x 10^ s, 24 s and 0.073 
s. Comparing this data with the Ti calculation of Ref. [3] shows that in spite of the strong suppression of relaxation in 
the Si QD due to a small deviation go— 9, the phonon- induced gi-factor fluctuation via excited states can control the 
phase relaxation (i.e., T2 < Ti) at T > 2 K. The slight anisotropy of the g factor [10] results in some T2 dependence 
on the magnetic field direction. However, this effect is expected to be small (less than 25%) in contrast to the strong 
magnetic anisotropy of T\ [3]. 

Qualitatively, another situation arises in the case of HFI modulation by the phonon- assisted transitions. The 
distinctive feature of this mechanism is an uncontrolled dispersion of the local nuclear field over the ensemble of 
QDs due to the random distribution of nuclear spins. This dispersion accounts for the fast (but partial) loss of 
initial electronic polarization of the QD aggregate without spin coherence loss [6,7]. Thus, we take into account 



the spin relaxation of a typical QD with a mean value of nuclear field dispersion (50„ = ay + l)Kni /Vqdi 
where / and n/ arc a nuclear spin and its concentration in a QD of volume Vqd] the dimensionlcss parameter 

K = Vqd J ^l^e (^1^ ^ \^g (Pfis equal 9/167r within the approximations of Ref. [3] and a is a constant in the 



HFI Hamiltonian Hhf = o-sl 



* [R 



for the nuclear spin / sited at a point R and an electron with its envelope 



function 'l'(r). If we set $7^ — fif = (557^ in Eq. (7), we readily find the following estimation 

T.-l = ^IiI^^)'>nF(l\. (14) 



This equation shows independence of T2^hf on the magnetic field in the case of HFI-induced transitions. An estimation 
of Eq. (14) (with an appropriate HFI averaging over ^^Ga, "^^Ga and "^^As [13]) for a GaAs QD with a typical size 
of 50 A x500 A x500 A and (5 = 2 mcV gives T2 = 34 s, 3.2 x lO'^ s and 1.0 x lO'^ s for T = 1 K, 2 K, and 4 
K, respectively. A similar estimation performed for a Si QD with the same dimensions and temperatures results in 
T2 = 4.1 X lO'' s, 3.8 X 10^ s and 1.2 s, respectively. One can see that for the considered cases of both GaAs and Si 
QDs, the contribution of HFI to phonon-assistcd relaxation is small in comparison to that of the g-factor modulation 
mechanism at i? = 1 T. However, the role of HFI modulation mechanism described in Eq. (14) should prevail at 
weaker magnetic fields {B < 0.1 T). 

The common property of Eqs. (12) and (14) are their inverse proportionality to the rate of phonon-induced tran- 
sitions between electronic states that qualitatively distinguishes this mechanism from those based on the spin-flip 
transition probabilities [4]. This also means that strong electron- phonon coupling may be favorable for slow decoher- 
ence (in a manner analogous to the D'yakonov-Perel process in the bulk materials). At the first glance, it would seem 
that can be enhanced limitlessly by increasing Tc- However, it will only be the case when the inequality Tc <C T2 
is not broken. An increasing Tc means that spin performs almost a full (or more) revolution during the procession in 
the fluctuating field between the sequential phonon- assisted electron transitions. Hence, in this a case we expect the 
dependence of Tj"^ according to the physical nature of the mechanism under consideration. 

Our formulas were obtained under the assumption of a two-level electronic structure (with regard to the orbital 
degree of freedom). In a similar manner, one can consider a many- level case. For example, the A^e excited electron 
levels close together with an identical g factors or HFI can be taken into account by replacing the tanh x term with 
a function (A^eC^ - e"^) / {N^e" + e"^) in Eq. (11). 

When considering the HFI as a mechanism of spin relaxation, we restricted our investigation to the SLR of single, 
but typical, QD with a mean configuration of nuclear spins. It is important to bear in mind that the relaxation 
measurements performed on a large ensemble of QDs correspond to averaging the relaxation curves over random 
T2jif instead of the single spin relaxation with a mean T2jif value. As a result, our mechanism would predict 
non-exponential decay of initial electron magnetization in a QD ensemble. 

In conclusion, we considered an eflacient mechanism of electron spin decoherence in a single QD due to the fluctuating 
precession of the longitudinal (with respect to direction of external magnetic field) effective magnetic field. Compared 
to the SLR controlled by the spin-flip transitions, our mechanism does not involve transitions between the Kramers 
doublets, which leads to eflfective spin relaxation characterized by a smooth dependence on the applied magnetic field. 
On the other hand, the mechanism under consideration reveal an exponential dependence on the temperature when 
T <C (5. This provides an opportunity to utilize the T2 temperature measurements as a tool for the study of QD's 
low energy structure. Although detailed estimation of the spin relaxation rate depends on the specific properties of 
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the respective QDs, the analysis on the typical cases illustrates an advantage of Si QDs for the quantum computing 
applications due to their very long relaxation time. 
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FIG. 1. Decoehrence factor F (5/2T) [Eq. (11)] as a function of temperature. 
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